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RELAXATION OF AN UNSTEADY TURBULENT BOUNDARY LAYER 
ON A FLAT PLATE IN AN EXPANSION TUBE 

By Roop N. Gupta and Robert L. Trimpi 
Langley Research Center 

SUMMARY 

An analysis is presented for the relaxation of a turbulent boundary layer on a semi- 
infinite flat plate after passage of a shock wave and a trailing driver gas-driven gas inter- 
face. The problem has special application to expansion- tube flows. The flow- governing 
equations have been transformed into the Crocco variables, and a time-similar solution 
is presented in terms of the dimensionless distance-time variable a and the dimension- 
less velocity variable /3. An eddy-viscosity model, similar to that of time-steady bound- 
ary layers, is applied to the inner and outer regions of the boundary layer. A turbulent 
Prandtl number equal to the molecular Prandtl number is used to relate the turbulent heat 
flux to the eddy viscosity. The numerical results, obtained by using the Gauss- Seidel line- 
relaxation method, indicate that a fully turbulent boundary layer relaxes faster to the final 
steady- state values of heat transfer and skin friction than a laminar boundary layer. The 
results also give a fairly good estimate of the local skin friction and heat transfer for near 
steady-flow conditions. 


INTRODUCTION 

When an unsteady flow occurs in the vicinity of a body, there is a time lag during 
which the flow characteristics adjust toward the new situation. This time lag results 
from different mechanisms in the external flow field and the boundary layer. The external 
flow field may have a quasi- steady character after a time required by the fluid particles 
to travel a few body lengths. This relaxation of the flow field is the result of wave- 
propagation mechanisms. The boundary-layer flow adjustment, however, is the result of 
molecular or turbulent diffusion, which is a much slower process than wave propagation. 

In an expansion tube (ref. 1) the interest lies in simulating the steady flow past a 
model. The distance-time plot of expansion- tube flow cycle is shown in figure 1. As 
pointed out earlier, the diffusion processes require time, yet the quasi-steady flow avail- 
able in the test facility exists for only a short time. It therefore becomes essential to 
ascertain that the establishment of the boundary layer (on the model) takes place in a 
time which is substantially less than the available test time. This situation is an asymp- 



totic one in which the departure from the final state is vanishingly small after a certain 
period of time. 

In references 2 and 3 results have been obtained for the establishment of a laminar 
boundary layer on a test model in an expansion tube. It is shown therein that an all- 
nitrogen time-dependent boundary layer on a flat plate relaxes to the Blasius value in 
about three-tenths of the distance (measured from the leading edge) traversed by the 
interface in the free stream. In practice, the boundary layer may not remain laminar. 
Under such circumstances it is of interest to know how the laminar results of refer- 
ences 2 and 3 are modified if the boundary layer is assumed fully turbulent from the 
beginning of its formation. In this study an analysis is presented for the wave-induced 
fully developed turbulent^compressible boundary layer (in region 5, fig. 1) over a semi- 
infinite flat plate in an expansion tube. The time and distance required to complete 
transition from laminar to turbulent flow are not considered herein. For quantitative 
evaluation of the boundary -layer properties, nitrogen has been used both in the test as 
well as in the acceleration gas region. The turbulent boundary layer is regarded as a 
composite layer consisting of inner and outer regions with separate eddy -viscosity 
expression for each region (in a manner similar to that employed for time-steady turbu- 
lent boundary- layer problems; see ref. 4, for example). The turbulent Prandtl number 
is assumed equal to the molecular Prandtl number and is used to relate the turbulent eddy 
conductivity to the eddy viscosity. The governing equations are treated in the Crocco 
transformed plane, and the conically self-similar solutions have been obtained by using 
the Gauss-Seidel line-relaxation method similar to the one employed in reference 2. 

SYMBOLS 

A ratio of shock velocity to free- stream velocity, U S /U Q 

a fraction of model length 

C = PP - 
(mP)w 

T-U/ 

c f local skin-friction coefficient, — 

i( P w u o 2 ) 

Cp specific heat at constant pressure 

c v specific heat at constant volume 
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G 


Blasius shear function, <^(0,(3) 


H dimensionless enthalpy difference 

H g Blasius enthalpy function, H(0,/3) 

Mirels enthalpy function, H(l,/3) 

H n dimensionless enthalpy difference behind shock in free stream, — 

o * n w 

h local enthalpy; also, step size in numerical method 

h Q local enthalpy behind shock in free stream 

Kj constant in eddy- viscosity model (eq. 40), 0.4 

K 2 constant in eddy-viscosity model (eq. 41), 0.0168 

k thermal conductivity 


krj* 

eddy conductivity 

L 

reference length 

M 

Mirels shear function, <p{l,0) 

M s 

shock Mach number 


free- stream Mach number 

N Pr 

C n^ 

molecular Prandtl number, — 

k 

Np r,t 

c p e 

turbulent Prandtl number, 

krp 

N st 

-q w 

Stanton number, 

"wVp T r ' T 
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n 


reciprocal ol velocity exponent in boundary layer, 


u 

Uo 



p pressure; also, iteration parameter 

q heat- transfer rate 

R gas constant 

R e x Reynolds number based on free-stream conditions, P e U 0 xyp e 

R w k Reynolds number based on L, p w U 0 L / /p w 

R w x Reynolds number based on x, P w U 0 x /^w 


R 


PooV 

Moo 


T temperature, K 

T q free-stream temperature in region between shock and interface 

T r recovery temperature 

T st stagnation temperature in free stream (having temperature T^) 

free-stream temperature in region between interface and expansion fan 

t time, measured so that at t = 0 shock wave is located at leading edge of 

plate 

t* time, measured so that at t* = 0 interface is located at leading edge of plate 

U e free-stream velocity in shock-fixed coordinate system 

U 0 free-stream velocity 

U g shock velocity 


4 



plate velocity in shock-fixed coordinate system 


velocity component parallel to plate surface in boundary layer 
velocity component normal to plate surface in boundary layer 


distance from leading edge along plate in x,t coordinate system fixed to 
plate 

distance from leading edge along plate in x*,t* coordinate system fixed to 
plate 

distance from foot of shock in shock-fixed x + ,y + coordinate system 
normal distance from plate 

location separating inner and outer regions of turbulent boundary layer 
normal distance from plate in shock-fixed x+,y + coordinate system 
conical coordinate in MR limit, x^Uot 
conical coordinate in BL limit, x* ju 0 t* 
dimensionless velocity, u/U 0 

Ue 

ratio of specific heats 

dimensionless time in MR limit, p w U 0 ^t/p w 




e 


o / 

dimensionless time in BL limit, P w U Q t / M w 

thickness of velocity boundary layer 

incompressible boundary- layer displacement thickness, 



dy 


eddy viscosity 


e = 1 +■£■ 


e = 


(V> 


3/5 


M N Pr,t k 


dimensionless y- coordinate 


^inc 


incompressible boundary -layer momentum thickness, 


dynamic viscosity 
kinematic viscosity 
density 
shear stress 
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shear function defined by equation (23a) 


dimensionless shear stress defined by equation (20a) 



Subscripts: 

e edge of boundary layer 
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I 

instantaneous value 

i 

inner region of turbulent boundary layer 

m 

count of steps in a-direction 

max 

maximum 

n 

count of steps in /3-direction 

0 

outer region of turbulent boundary layer 

t 

total condition 

w 

evaluated at wall 

oo 

free stream 

Superscripts: 

t 

fluctuating component 

+ 

quantities in shock-fixed coordinate system 

Abbreviations: 

BL 

Blasius 

MR 

Mirels 

N 2 

nitrogen 


BACKGROUND 

There does not appear to be any existing literature pertinent to the expansion-tube 
turbulent boundary-layer problem considered herein. For a plate located ahead of the 
secondary diaphragm, the problem is analogous to the shock- tube problem analyzed in 
references 5 and 6. For a plate located behind the secondary diaphragm, however, an 
accelerating-gas boundary layer is first established on the flat plate. This boundary 
layer is subsequently swept over by the interface separating the test gas and the acceler 



ating gas at two different temperatures. Finally, a test-gas boundary layer is established 
close to the leading edge of the plate. Whereas Mirels* problem (ref. 5) is steady in a 
shock-fixed coordinate system, the present problem is essentially time dependent. Thus, 
no semiempirical formulation (as used by Mirels and others) can be made. 

MATHEMATICAL FORMULATION 

Turbulent Compressible Boundary- Layer Equations 

The basic differential equations describing the two-dimensional flow (with instan- 
taneous values of velocity, density, and enthalpy) in a boundary layer on a flat plate (with 
the assumption of pj = Constant) are (ref. 7): 

Continuity: 


S P T 9(pu) T a(pv) T 

— + ± + ± = 0 

8t 9x 9y 


Momentum: 


1 9u T 8 u t 9u- 


Pj — + Uj — 3 + Vj — 5 = — (T yx ) 

I \9t 1 0x 1 ay/ 9yW X /l 


Energy: 


where 


/ 9h T 8h T 9h T \ n . 9 u t 

Pt — 7 + u i — - + Vj = — (qJ) + (r ^ 

I \9t 1 9x 1 By 9yW/i V^/iay 


( T y x )j = ( T xy) ! " M Ity 


and 


(q\ * 

\ y h Np r 9y 


( 1 ) 


( 2 ) 


(3) 


( 4 ) 
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Following Van Driest (ref. 7) the instantaneous values are separated into average 
and fluctuating components: 


Uj = u + u' 

P 1 = P + P' 
hj. = h + h’ 

Tj = T + T' 

Vj = v + v’ I 



T + T' 

yx yx 


( 5 ) 


( q y)i q y + q y 

) 


The following averaging procedure is carried out with respect to time at a fixed 
point in space for a time-dependent mean flow. If At is a time interval much greater 
than the characteristic time scale associated with the fluctuations in a quantity f(x,y,t), 
then the average value of f(x,y,t) at a given point (see fig. 2(a)) is described by the 
equation 



( 6 ) 


In case f (x,y,t) does not depend on the averaging interval At and physical time t, the 
mean value is said to be independent of time (see fig. 2(b)). 

Introducing the expressions defined by equations (5) into equations (1) to (3), aver- 
aging with respect to time as defined by equation (6), and using the thin-boundary-layer 
assumptions, the following equations are obtained (see ref. 7 for details): 

Continuity: 

9£ + 3(PU) + 8(pv) _ o (7) 

9t 9x 9y 
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Momentum: 


0U 3u ~ 9u 9 [, ; 9u\ 

p — + PU — + pv — = — pe — 

H at 3x ay ay \ ay j 


(8) 


Energy: 


ah 9h - 3h a / fi * 9h\ , -W 

at ax ay a ylNp ay J y>y t 


0) 


where 


v = v + 


pV 


( 10 ) 


6 = 1 + — 
p 


p Np r>t k 


( 11 ) 

( 12 ) 


The eddy viscosity has been defined as 


e = 



(13) 


and the eddy conductivity as 


lOp = -P 


v'h' 

aT/a y 


(14) 


The static turbulent Prandtl number is defined as 


N 


Pr,t 



(15) 
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The following perfect-gas law is also introduced: 


P = 



( 16 ) 


Other assumptions made are similar to those made in reference 2; namely, 

(1) Since no general time- similar solution is possible, the problem is treated in two 

conically similar limits: (a) when the time-lag between arrival of the shock and the inter 

face at the leading edge of the plate is very large (Blasius or BL limit), as shown in fig- 
ure 3; (b) when this time lag is negligible (Mirels or MR limit); this limit is depicted in 
figure 4. 

(2) The flow external to the boundary layer is not affected by the boundary-layer 

flow. 

(3) The interface in the inviscid flow is considered to be thin. 

(4) The wall temperature is assumed to be constant. 

Transformed Plane 

As pointed out in reference 2, the Crocco system of variables is suitable for treat- 
ing time-unsteady problems such as the one to be treated here. In this system the inde- 
pendent variables are 


x = x 
u e u(x,y,t) 
t = t 


and the dependent variables are 


9u 




T = (pi + €> 9 ^ = T (x,u,t) 




h = h(x,u,t) 


(17) 


(18) 
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In the conical coordinate system the independent variables (eqs. (17)) can be written as 


a = 


V 


j3 = 


= U 0 / 


(19) 


- W* 


Y = 


W_ 

M 


w 


and the dimensionless dependent variables as 


<P(a,$,Y) = 

P w U o 


(20a) 


H(a,j3,y) = 


h - 




(20b) 


'w 


With the introduction of equations (17) to (20), equations (8) and (9) become, respectively, 


3 2t 






y Bq! By 


a - j3 9(Ce) 8(Ce) 

y 8(2 9y 


( 21 ) 


0 2 


U Q t 8 [ 6 1 9H 1 

0 i3\£ Npr 0 0 


1_£ . At 8 0 8 H _ C J 8H _ a - g 
^Np r 6 j 9)3 9)3 \^y y 


9H’ 

9a) 


( 22 ) 


In obtaining equations (21) and (22), the continuity equation, equation (7), has been utilized. 

To obtain self-similar solutions in the Crocco system, the following relations are 
assumed for <f>, H, and C: 


0(a,j3,y) 


<«r> 1/5 


(23a) 
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thus H = H(ai,0) and it follows that 


0C 

8 y 


C = C(a ,$) J 


(23 c) 


If the thermal- wall conditions are time dependent, then the assumptions contained 
in equations (23) cannot be made (see ref. 8 for details). In that case, the presence of 
three independent variables (o!,jS,y) instead of two will substantially complicate the compu- 
tation of the solutions. Presumably a step-by-step forward integration in the positive y 
direction may.be used to treat such problems. 

Equations (21) and (22), with the use of equations (23), can now be written as 


Ce 


^ («y ) 3/5 


. (h 

a{a - |8) — + j3 ~ 
8<* 5 


<P 


, -,3/5 

(ay) 1 


«(a - P) 

da 


I 

9 y 


(24) 




U o 2 1 8 /e 1 &h\ 

hw 9/3\e Np r 90 J 



90 9/3 


Ce 



a(0 - a) — 
9a 


(25) 


Further, in order to make equations (24) and (25) self-similar, it is necessary to 
assume the following relation for e : 

€(ay,a,0) = (o-y) 3 ' /5 €(Q!,0) (26) 

The validity of this equation will be discussed later. 

The Incompressible data which are available for the core region of pipe flow and the 
outer regions of boundary layers (ref. 9) indicate that Np r ^ ranges between 0.7 and 0.9. 

For compressible flows, where very little data are available, Meier and Rotta (ref. 10) 



found for 1.75 = = 4.5 that Np r ^ 


increased above unity for y* < 50 


where 


and ranged between 0.8 and 0.85 as the outer edge of the boundary layer is 

approached. 

In this analysis the molecular Prandtl number Np r and the turbulent Prandtl num- 
ber N Pr t will be assumed, for simplicity, to have the same constant value of 0.8. For 
this case 



€ = € (27) 

In view of the relations in equations (26) and (27), equations (24) and (25) finally yield 


<p 2 — + - Cefi<f> = a(j3 - ot)e (c | - - <p - Qf((3 - a)C<p ~ 
V 0 ^2 5 Y \ da da J da 


(28) 




l w 


N Pr 9 / 3 2 


IN 


Pr 


9/3 9/3 


a ( 0 - a) 


9H 

da 


(29) 


BOUNDARY CONDITIONS 

The various regions used in the specification of boundary conditions have been 
defined in reference 2. The same definitions will be used in this section. In the trans- 
formed plane the regions of interest in the BL (or Blasius) limit and the MR (or Mirels) 
limit have been shown in figures 5 and 6, respectively. 

Shock Region (M) and Limited Interaction (LI) Region (Contained 
Between a* = 1 and a = 1 in Figs. 3(b) and 5) 

Boundary Conditions 

First the shock region is considered. The separation- of- variables approach used 
for the laminar case (ref. 2) does not work for the turbulent case. In order to overcome 
this difficulty, the solutions for the shock region have been obtained in the shock- fixed 
coordinate system (fig. 7). The shock-region boundary layer in this coordinate system is 
time steady, and the solutions can easily be obtained. The following relations exist 
between the plate-fixed and shock-fixed systems: 
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(30) 


x + x+ = U~t 


Uo + U e = U s 


<t>(a,i 3) = al/5 (4 . - -. 1 ).- f M + (j3) 
(A - *) 1/5 


(31) 

(32) 


H(a,/3) = H^(/3) 


(33) 


The equations governing M + (0) and H^(/3) are (see appendix A for details): 


M + + 4 £ c -+ = 0 


d|3' 


,2 5 


(34a) 


M 


1 o 2h m 


N Pr dj 3 2 


, ( U e - M 


A w 


lN Pr I d/3 d/3 


(34b) 


The boundary conditions on M + (/3) and H^(/3) are 




dNT 


d/3 


M 


= 0 


J3=0 


0=1 


= 0 


(35a) 


H 


M 


= 0 


0=0 


H 


M 


= H„ = 


h o " 


0=1 


0 h 


w 


(35b) 
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The solutions of equations (34a) and (34b) are used in relations (32) and (33) to pro- 
vide the boundary conditions on <p and H at a = 1 for the I region in the MR limit. 
The I- region boundary conditions on (p{a * = 1,0) and H(a* = 1,0) for the BL limit are 
obtained by solving the steady-state forms of equations (28) and (29) (see eqs. (36a) 

(36b)) subject to the boundary conditions (37a) and (35b), respectively. 


Leading-Edge Boundary Conditions 

The boundary conditions at a* = 0 (BL limit) and at a = 0 (MR limit) are pro- 
vided by the solution of the steady turbulent boundary-layer equations for a flat plate. At 

B<P\ 


a = 0, <p = G(0) and, for an asymptotic approach to the steady state — j - 

da 

momentum equation (eq. 28) for such a case yields 


= 0. The 


a=0 


+1 Ce0 = 0 
d 0 2 5 


(36a) 


Similarly, by taking H| n = H r (0) and — 

QE=U Ci dcx 

simplifies to 


= 0 the energy equation (eq. (29)) 


a = 0 





Pr 


.itefEs.o 

I d/3 d/3 


(36b) 


The boundary conditions on G(/3) and H G (/3) are 


A 


*£ =o 

d/3 

P > 13 = 0 


G U = ° 




(37a) 


H, 


H, 


0=0 


= 0 


0=1 


= 0 


(37b) 


The BL- limit boundary conditions at a* = 0 can be obtained by replacing a by 
at* in assumptions leading to equations (36) and (37). 
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Wall Boundary Conditions 

For zero slip, impermeable wall, and constant wall temperature, the MR-limit 
boundary conditions on cj> and H are 

££(a,/9 = 0) = 0 

d/3 

H(a,/3 = 0) = 0 

The wall boundary conditions for the BL limit may be obtained from equations (38) by 
replacing a by or*. 

Boundary- Layer Edge or Free-Stream Condition 
For MR limit the free- stream boundary conditions on (p and H are 


^0 = a = lj 


( 38 ) 


<Ma,/3 = 1) = 0 

H (a,/3 = 1) = 0 




(39) 


Once again, the BL- limit boundary conditions may be obtained by using a* in place of 
a in the equations (39). 

In equations (39) the free-stream temperature (between the interface and the expan- 
sion fan) has been taken equal to the plate temperature. This assumption is appropriate 
for expansion- tube flow duplication studies (ref. 11). For simulation studies equa- 
tions (39) should be modified for relevant values of H in free stream. 

The values of the parameters A, U Q 2 /» w , and T o/ T w required in obtaining the 
boundary conditions and for solving the governing equations have been obtained from the 
inviscid shock-tube relations (see appendix B). 


TURBULENT TRANSPORT MODEL 

Following references 4 and 9, the turbulent boundary layer is treated as a composite 
layer characterized by inner and outer regions. In the inner region, an eddy viscosity 
based on Prandtl's mixing- length theory is used; in the outer region, a nearly constant 
eddy viscosity based on the Clauser model (ref. 9) is used. It is exactly constant when 
the flow is incompressible and without heat transfer. The intermittent character of the 
outer region is accounted for by using an intermittency factor. 
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Viscosity in the Inner Region 

The eddy viscosity for the inner region referenced to the molecular viscosity is 
obtained from (ref. 4) 




where A = 26v J — and v is the mean value of v taken over the viscous sub- 


layer. As an approximation the ratio of is assumed to be unity. 

Viscosity in the Outer Region 

The ratio of the eddy viscosity to the molecular viscosity in the outer region can be 
expressed as (ref. 4) 


e \ rr o.* - 

J 2 M 5 inc y 


where 6. is the incompressible displacement thickness 




and the intermittency factor y is given by 


1 - erf 5/i - 0. 78 


In this expression 6 is the thickness of the boundary layer. 


Matching Procedure 

The criterion of the continuity of eddy viscosity gives the following relations for the 
matching at the boundary between the inner and outer regions: 
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pK^y 2 




*Vo * 

u °inc 


y 


where y m is determined from 


(° = y = y m ) 


(y > y m ) 

J 


( 44 ) 



(45) 


Transformed Models 

From equations (18), (20), and (23) the following expression may be obtained 


3u _ Pw U o <k 

“ + £ («^ 1/6 


which, upon integration, gives 


y = 


p w u 0 -'o 4> 


Equation (47) may finally be written as 


(46) 


(47) 


y = 




(48) 


Substitution of expressions (46) and (48) in equation (40) results in 



K^I 2 
C 2 (H + l) 3 






0 


(49) 
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where the following relations have also been used (from perfect- gas relations and 
boundary -layer approximations): 



Equation (49) is quadratic in (e/p)^ and, therefore, has two solutions for (e/p^. 
It is easy to see, however, that (e/ pj^ goes to zero as y goes to zero only if the posi- 
tive root is taken. Hence, the relation for (e/p^ is obtained as 



Similarly the eddy-viscosity equation (eq. (41)) for the outer region can be written 
in the transformed plane as 


=1 — 

^'o Cw 


— ~ 2 ^ _ ( Rw > x ) 1/5 1 " erf 5 ( r ~ 

(H + l) 2 2 ' ' 1*0.99 


- 0.78 
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where 



( 53 ) 


Now, when the exponential term of equation (40) is expanded and a logarithmic 
approximation is used for the velocity profile, it is found that (e/V)j varies as y 3 . Of 
course, such a logarithmic approximation is assumed to be valid only for small values of 
y lying in the region between the laminar sublayer and the region where the classical 
logarithmic law with a different slope is often applied. Hence, from equation (47) it may 
be said that (e/g^ varies as (R^x)^ 5 - 

In order to make equations (24) and (25) self-similar, relation (26) was employed. 
The relation in equation (26) gives the following Reynolds number dependence for ( e//i )^: 




(54a) 


and the further approximation is introduced 


£i(a,j3) 




(54b) 


where (e/ju)^ is evaluated from equation (51) and R^l is evaluated at x = L, which is 
the average length of the model. This approximation is not valid for x « L but is used 
since the basic assumption of turbulent flow is on no firmer grounds. 

In order to make equation (52) compatible with relation (26), equation (52) is approx- 
imated as 


so that 





~ 1 


*2 


C , 2 , .2/5 2 

(H + 1) (R W)L ) 


1 - erf 5/ 


l 0.99 



(55) 


( 56 ) 
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Before closing this section it should be pointed out that the assumptions made in 
obtaining expressions (54) and (56) are further additions to the original assumptions con- 
tained in expressions (40) and (41), which themselves are also approximate. For lack of 
anything better, however, these relations are employed and their validity checked against 
some time-steady experimental data. The validity of these results for the unsteady flow 
still needs experimental substantiation of equations (54) and (56). 

SOLUTION OF THE GOVERNING EQUATIONS 
Numerical Procedure 

The numerical solution of the boundary -value problem (eqs. (28) and (29)) formu- 
lated in this work is sought by using the Gauss-Seidel line-relaxation method. The finite 
differencing scheme and the method of linearization used here closely parallel that of 
reference 2. 

By using the central difference approximations (see ref. 2) equations (28) and (29) 
may be written as follows for solutions along lines of constant 0 (the superscript p 
refers to the iteration number): 


^ m ( a m ■ 
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2 
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In equations (57) and (58) the nonlinearities have been removed by using values of 
the dependent variables from the previous iteration in the coefficients. Each of equa- 
tions (57) and (58) forms a set of -lj simultaneous difference equations for a given 

value of n when m varies from 2 to 1/h. Since the coefficient matrix of each of 
these equations is of the tridiagonal form, the Thomas algorithm of reference 12 may be 
employed. 

The numerical solution is started by calculating the fluid properties and establish- 
ing the inner and outer regions from laminar <p and H distributions. The momentum 
equation (eq. (57)) is now solved iteratively for a new <p distribution along lines of con- 
stant /3 for the whole flow field until 



* _ CL 

where is of the order of 10 . This value of p is taken as P max - 

Next, equation (58) is solved iteratively in the same manner P max times. The 

P mnv P m a -y 

improved <p and H distributions are now used to calculate the new fluid 

properties and establish inner and outer regions. With these values the momentum and 
the energy equations are solved again. This process is continued until convergence to 
the desired accuracy 


<ft J (q,0) - <ft J %b0) 

is obtained. Here e Q * is a small number of the order of 10” 5 and J is the iteration 
number. 

The convergence rate of the solutions can be improved substantially if equations (57) 
and (58) are solved first for a coarse mesh (h = 0.1), and then with linear interpolation 
these coarse solutions are used as initial guesses in the iteration with fine mesh 
(h = 0.01). A computer flow chart for the procedure used is given in figure 8. 
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Calculated Boundary- Layer Parameters 

Once the <p and H distributions from the numerical solutions are obtained, the 
boundary -layer parameters, such as local skin-friction coefficient, Stanton numbers, and 
boundary -layer thickness, can easily be obtained. Some of these parameters are given 
by the following equations. The local skin-friction coefficient is given by 


C f = I 


w 


f w 


(59) 


where 


T w ' ^w 




From equations (20) and (23) 


T w _ ^w 

p w u Q 2 


(60) 


Hence, 


c f = 


2 _K 

(«r ) 1/5 


or, 



The heat transfer at the wall can be obtained from 
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In terms of the transformed variables, this result can be further written as 






N Pr 


Defining the Stanton number as 


N St = ' 


4 W 


p w U o c p( T r ' T w) 


( 63 ) 


(64) 


equation (63) becomes 


N = 1 / 8H\ 1 

St N Pr (a^) 1 / 5 v^/w T r 1 


(65) 


L w 


where the recovery temperature T r is obtained from the ratios 


r = 


Tj. Too 
T st " T °° 


- (N Pr ) 


1/3 


( 66 ) 


1/3 

The assumption of r = (Np r ) has no firm basis for such unsteady flows and is 
probably inaccurate. It is used here simply to permit expression of the results in terms 
of a Stanton number. With 


Tst - I'oo + 


U„ 


2c r 


and 


Too — T w 
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expressions (66) yield 




1/3 V 
2h w 


(67) 


From equations (65) and (67) the following expression is obtained for the Stanton number 



2h w ‘K fm\ 
U o 2 (Np r ) 4/3 ^' 5 ^ 


(68) 


Inversion to Physical Plane 
From the definition of the transformed variables 

r = p U 2 £ — 

MV o 1 /c 

(<*5} 1/5 


or 


so that 


/ 4 9 \ 1/5 

9u _ | P W U 0 P W I 0 

a y \ x / m + e 




-* 1/5 

.1/5 ^0 


( p w U o ^Wj 


Now 


(69) 


(70) 


p = Cp w p wf-^ (H + 1) 


(71) 
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Hence, from equation (70) the dimensionless ordinate rj may be written as 


rj(a,/3) = 


4 4 

x Pw ^w 

U„ 


1/5 


Jq cp 


(72) 


y - 1 h w 

r p n 


where e is related to e by the expression in equation (26). 

By evaluating the integral appearing in equation (72) at fixed a for various £ 

values, the profiles of r}(a,r^-^j at various a values can be obtained. These may be 

used to obtain velocity profiles — ( 77 ) for various values of a. 

Uo 

The temperature in the boundary layer may be obtained from 


h 

T = — (H + 1) 
C P 


(73) 


From equation (72) 


(* C(H + Dc ^ 

I = -J— , t (74) 

6 ^0.99 r°- 99 C(H + l)g 

J 0 (p 


may also be obtained. 


DISCUSSION OF RESULTS 

Evaluation of Some Numerical Results and Comparison With Experimental Data 

Some of the numerical results for skin-friction drag for a fully turbulent boundary 
layer on a flat plate are compared herein with the experimental data of Moore and 
Harkness (ref. 13). The experimental results were obtained on a sharp- leading-edge flat- 
plate model with a floating-element-type balance for 1 x 10^ < Roo )X < 2 x 10®. The test 
conditions were as follows: 
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M to = 2.8 


Pt oc = 0-9997 X 10 6 N/m 2 
T t0O = 311.1K 

L J 

T 

— — = 0.947 
T t~ 

The governing equations, equations (28) and (29), reduce to the time-steady flat- 
plate fully turbulent boundary -layer equations for a = 0. The numerical results of these 
equations, solved to obtain the local skin-friction coefficient, are compared with the 
experimental data in figure 9. 

The agreement is good at high Reynolds numbers and fair at lower Reynolds num- 
bers. Note that the value of the reference Reynolds number, R T , used in the compu- 

7 w,1j ^ 

tations is 1 x 10 . For the range of R e x considered in figure 9, x/L varies from 0.4 
to 4. 


Time-Dependent Results 

The skin-friction, heat-transfer, velocity, and temperature profiles have been 
obtained for a wave-induced compressible turbulent boundary layer on a flat plate after 
passage of an interface with application to expansion-tube flows. The gas considered on 
the two sides of the interface is ’’perfect” nitrogen. The results presented correspond to 
a free-stream velocity of 6 km/sec with a shock Mach number of 20.8. A reference 
Reynolds number of 3.78 x 10® per meter is employed in expressions (54) and (56). This 
value corresponds to the following conditions: 

p w = 10 3 N/m 2 

U Q = 6 km/sec 

T w = 300 K 

M w - 1.785 x 10 5 N-sec/m 2 

The molecular and the turbulent Prandtl numbers are assumed to have the same 
constant value of 0.8. The plate as well as the free stream behind the interface is 
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assumed to be at the same temperature, 300 K, for this analysis; this assumption is per- 
tinent to flow- duplication studies (ref. 11). 

Figures 10 and 11 give the local skin-friction coefficient for the MR and BL limits, 
respectively. In both of these figures the skin-friction coefficient relaxes to within 6 per- 
cent of the steady-state value at a = 0.6. Shown in the following table are the wall skin- 
friction coefficients computed by the Spalding-Chi correlation (ref. 14), by the Van Driest 
method (ref. 7), and by the Schultz -Grunow method (ref. 15) employing the Eckert refer- 
ence temperature approach to account for compressibility. It may be noted here that the 

present problem produces c^R w?x ^ 1//5 = Constant for the steady state, whereas the other 
methods do not. 


For a = 0 (steady-state results) 

^e,x 

Spalding-Chi 
(ref. 14) 

Schultz - Grunow 
(ref. 15) 

Van Driest 
(ref. 7) 

Present 

c f( R e,x] 

/r T,„ 

' for — = 1; M e ~ 17; T w - 300 K 
T e 

1.17 x 10 6 
1.17 x 10 7 
1.78 x 10 8 

6.054 x 10 

6.880 

8.420 

-3 

11.26 x 10' 3 

10.40 

10.67 

10.55 x 10" 3 

10.80 

12.00 

J 

> 9.22 x 10" 3 

N St( R e,x 

) 1/5 for ^ = 1; M e ~ 17; T w = 300 K 

1.17 x 10 6 
1.17 x 10 7 
1.78 x 10 8 

3.27 x 10" 

3.67 

4.44 


6.43 x 10“ 3 

5.75 

5.76 

1 

5.86 x 10" 3 

5.86 

6.43 

J 

> 4.8 x 10 3 


The Stanton number distribution for the MR and BL limits is also shown in fig- 
ures 10 and 11, respectively. These figures indicate that the wall heat transfer relaxes 

to (within 6 percent) its final steady-state value at a - 0.5. The high (|Sj values in 

\ P/ w 

the vicinity of a - 1 result from the imposed temperature boundary condition. Also 
shown in the preceding table are the Stanton number values computed by using the pre- 
viously quoted skin-friction relations from references 7, 14, and 15 and the Reynolds 
analogy. 

Since this work was started with the aim of determining how a wave-induced bound- 
ary layer would relax if it were assumed turbulent instead of laminar from the beginning 
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of its formation, a brief comparison with the laminar results is of interest. It was shown 
in references 2 and 3 that for a laminar boundary layer the Stanton number and the local 
skin-friction coefficient relax to within about 6 percent of the steady-state value of 
a. - 0.3 (for the same free-stream velocity and temperature values). The present work 
indicates that a fully turbulent boundary layer relaxes faster to the time-steady values as 
compared to a laminar boundary layer. This result should be expected from the enhanced 
transfer of momentum and energy in a turbulent boundary layer. The result is also sup- 
ported by the work of Trimpi and Cohen (ref. 6). They showed that the slope of the 
boundary-layer characteristic which bounds the steady-state is controlled by the ratio 

6* fa. . The Blasius value of this ratio for the laminar case (n = 1) is 2.6 and for the 
inc/ inc 

turbulent case (n = 7) is 1.29. These values imply that the slope of the boundary- layer 
characteristics is smaller when compared to the laminar boundary layer. Hence the tur- 
bulent boundary layer should relax faster as indicated by limiting values of a *0.39 
and 0.77 (reciprocals of 2.6 and 1.29) for the laminar and fully turbulent cases, 
respectively. 

Figures 12 and 13 show the velocity and temperature distributions through the 
boundary layer for the MR limit. Note that at a = 0.6 although the velocity profile has 
nearly reached its asymptotic shape throughout the entire boundary layer, the outer part 
of the temperature profile is still in the process of significant relaxation. 

The wave-induced boundary- layer behavior can be interpreted in two different ways 
under the conical assumption. The present results are presented for the development of 
the boundary layer with time at a given x position on the plate ("fixed x"). A descrip- 
tion of the boundary layer at any point in time ("fixed time”) is often desirable. The 
Crocco conical coordinates, (employed here) readily permit interpretation of results for 
either point of view. For the fixed-time point of view care should be taken f in obtaining 
the results at the leading edge and at the shock wave since the boundary- layer equations 
are singular at both of these locations. When the fixed x point of view is taken, the 
product ay produces a Reynolds number which permits interpretation for fixed values 
of x, and in such a case a becomes a time variable. For the fixed-time point of view, 
Y is fixed and a becomes a position variable. 

CONCLUDING REMARKS 

Crocco variables have been employed to obtain time-similar (or conical) solutions 
of the wave-induced flow over a flat plate in an expansion tube. The boundary- layer flow 
was assumed to be fully turbulent. Results presented correspond to two limiting cases in 
which the flow external to the boundary layer can be considered similar in time or 
"conical." All intermediate cases should lie within the range of these two limits. 
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From these results, it is seen that the fully turbulent boundary layer relaxes faster 
to the time- steady state when compared to the laminar boundary layer. For the cases 
analyzed herein the turbulent boundary layer relaxes to the final steady-state values 
(within about 6 percent) of skin-friction and Stanton number over a length, measured from 
the leading edge of the plate, equal to about one-half of the distance traversed by the inter- 
face (separating cold and hot sections of the flow field) in the free stream. 

The assumptions and approach were checked against steady-flow data and found to 
give a fairly good estimate of the local skin friction and heat transfer. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 16, 1974. 



APPENDIX A 


SOLUTION OF THE SHOCK- REGION GOVERNING EQUATIONS 

Since the separation-of-variables technique, used in the laminar case (ref. 2), was 
apparently not successful in the plate-fixed coordinate system, the solution for the shock- 
region (M) governing equations has been sought in the shock-fixed coordinate system. 
These results were later transferred to the plate-fixed system. (A similar treatment for 
the laminar case is given in ref. 16.) 

As depicted in figure 7(a), the shock wave is moving with a velocity U g and the 
fluid behind it with a velocity U Q . If the shock wave is fixed by subtracting U s from 
every part of the flow field (including the wall), the flow shown in figure 7(b) results. 

This transformation to a shock- fixed coordinate system removes the unsteady character 
of the problem (for the portion of the flow field which is not influenced by the leading edge) 
present in the plate-fixed reference frame. The problem can now be treated by analyzing 
the steady- state turbulent boundary -layer equations for a flat plate. 

Continuity: 


8(P + U^ , 8(p + v+) _ 0 
0X + a y + 


with 


+ pV 
v = v + t — L - 


(Al) 


Momentum: 


<p v> ^ = -yW + ' 

0x + 8y + 0y + \ 0y + y 


(A2) 


with 


i + = 1+1- 
M 
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Energy: 


(p + u + )^ + (p + V + )^ = € + + —V i a+e “ + (H) 

8x + 0y + 9y + \^Pr 9y + / \dy + J 


(A3) 


with 


.+ n 


M 1N Pr,t 


The following definitions of the dependent and the independent variables are now 
employed in the Crocco system: 

Dependent variables: 


T+ =■ U* * e + )^ 

\ >9y + 

h + = h + [x + ,u + ) 

Independent variables: 




(A4) 




x + = x + 


{i + ^ u + (x + ,y + ) 




In terms of the Crocco variables, equations (A2) and (A3) may be written as 


(A5) 


9 2 t + _ + 8 + 

o +u — + P Ji -T- =0 

- + 2 8x V T / 


(A6) 


(p + G + )M + e + ~= r + 
8x + 


1 + 


a r l 9h + ' 

0u + \e + ^Pr 9u + 


1/ 1 r -UT+ 9h + 
+ 2\N r r+ “ J Q .-.+ *r.+ 


'Pr e 


Su 8u 


(A7) 
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The following dimensionless variables are now introduced: 







(A8) 


The momentum equation (eq. (A6)) and the energy equation (eq. (A7)) then become, 
respectively, 


[1 

V u c 


i — —El + - t ^- —fiL i 

U„A 2 9/3 2 Pw U e d x + \6 + ' 


= 0 


(A9) 


juVe + ~/ Tr TT \ 2 9H. T +* 

^T 7 T e 'N — -0 


.+ 2 tt 3 
p w U e 


( u e - Uw)‘ 


1 BH' 


l w 


8 P e + N Pr 9 B 


— + 


1_ _ A + 9^_ 9H 


iNpr g + 


9/3 90 


(A10) 


where 


t.ZL.XZL' 


u. 


u c 


The following forms for <p + and e + are now introduced so that the momentum 
equation (eq. (A9)) becomes independent of x + : 

A/5 

M + (0) (All) 


0 + = fi - M 

\ UJL+ 


e /\p;u e x^ 
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i + = 


3/5 

fW ?+ 






w 


For the energy equation (eq. (AIO)) it is assumed that 


(A 12) 


H = H(j3) (A13) 

By using expressions (All) to (A13) in equations (A9) and (AIO) the following form of the 
shock-region governing equations is obtained: 


M + + | j3Ce" + - 0 

d/3 2 5 


M 


(u e ■ u w y 


w 


1 d 2 H | 
N Pr d/3 2 


i N 


Pr 


. i ) dH = o 

r d/3 d/3 


(A 14) 


(A 15) 


In obtaining equation (A15) it is assumed that Np r ~ Np r = 0.8. The preceding assump- 
tion implies e + = e + . 

It should be noted here that the nondimensional local enthalpy H has the same 
value in the plate and shock- fixed reference frames; however, the shear function M + (/3) 
is to be related to the shear function <p{x,i 3) in the plate-fixed reference system. 

Now, <p by definition is 


<p = 


PwUo 



or 


(A 16) 



(A 17) 
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In equation (A17) the minus sign has been dropped. It is also noted that 0 is of 
opposite sign to 0. 

The following relations exist between the plate-fixed and the shock-fixed reference 
frames: 



(A18) 


From equations (All) and (A17) the following relation is obtained for shock regions 
after substituting for x + and U e : 


(A19) 

therefore , 

r ii 1/5 

0(or,j3 ) = (ay) 1/5 0(or,fty) = M + (ffi (A20) 

A - a 

For 0 = 1 equation (A20) gives 

0( 1,0) = (A - 1) 3/5 M^) (A21) 
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APPENDIX B 


SHOCK PARAMETERS IN AN IDEAL GAS 


To be consistent with the rest of the analysis and for the sake of simplicity, the fol 
lowing properties of the gas are assumed: 

(1) The ideal equation of state. 

(2) Cp = Constant, c v = Constant. 

With these assumptions, the following relations are obtained from reference 17: 


A 


U 

U 


s _ y + 

r 2 


i 



h ° (y + l) 2 T 0 \ 



(Bl) 


(B2) 



- 1 ) 

(y + l) 2 




y) 


(B3) 


where Tj is the temperature of the gas in front of the shock. The temperature Tj 
has been taken as T w (=300 K) in this analysis. Further, 



h w 1>Q h^y hQ T j 


(B4) 


since c p is assumed to be constant and T^ = T w . 
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Figure 1.- Distance-time plot of expansion-tube flow. (Taken from ref. 1; 
circled numbers are explained therein.) 



(a) Time-dependent mean flow. (b) Time-steady mean flow. 

Figure 2.- Schematic diagram showing time-averaging method employed. 
(At, averaging interval; t c , characteristic fluctuation period.) 
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(b) Schematic representation of flow field in MR limit (physical plane). Time delay 
represented by shaded region is assumed to be vanishingly small. 


Figure 4.- Schematic representation of MR limit. 
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(a) Schematic representation of flow field in a*, $-plane. 
Figure 5.- Blasius (or BL) limit. 
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(b) Specification of boundary conditions on <p in I region. 


Figure 5.- Continued. 





►f* 

CT> 


H 



(c) Specification of boundary conditions on H in I region. 
Figure 5.- Concluded. 
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(a) Schematic representation of flow field in a, plane. 
Figure 6.- Mirels (or . MR) limit. 
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(b) Specification of boundary 
Figure 6.- 



on (p ■ in I region. 




(c) Specification of boundary condition on H in I region. 
Figure 6.- Concluded. 
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y 



u w = o 

(a) Plate-fixed coordinate system. 



(b) Shock-fixed coordinate system, 
'igure 7.- Shock- region boundary layers. 
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Figure 8.- Flow-diagram for solving turbulent boundary -layer equations. 
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Local skin friction coefficient 



Figure 9.- Variation of local skin-friction coefficient with Reynolds number. 






Figure 11.- Distributions of local skin-friction coefficient, Stanton number, 
and gradient of dimensionless local enthalpy; BL limit. 
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Figure 13,- Temperature distribution through boundary layer for 
M s = 20,8 in nitrogen; MR limit. 
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